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Screening of a macroion by multivalent ions: A new boundary condition for the 
Poisson-Boltzmann equation and charge inversion. 

V. I. PereH and B. I. Shklovskii 
Theoretical Physics Institute, University of Minnesota, 116 Church St. Southeast, Minneapolis, Minnesota 55455 

Screening of a macroion by multivalent counterions is considered. It is shown that ions form 
strongly correlated liquid at the macroion surface. Cohesive energy of this liquid leads to strong 
additional attraction of counterions to the surface. Away from the surface this attraction is taken 
into account by a new boundary condition for the Poisson-Boltzmann equation. This equation is 
solved with the new boundary condition for a charged flat surface and a long cylinder. For a cylinder 
Onsager-Manning theory looses its universality so that apparent charge of the cylinder is smaller 
than their theory predicts and depends on its bare charge. It can also vanish or change sign. 

PACS numbers: 77.84.Jd, 61.20.Qg, 61.25Hq 
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I. INTRODUCTION 

Many polymers are strongly charged in a water solu- 
tion. Lipid membranes, DNA and actin are the most 
famous examples of such biological polyelectrolytes. In 
a water solution polyelectrolytes are screened by smaller 
ions of both signs. Correct description of the screening 
of polyelectrolytes is tremendously important for calcu- 
lation of properties of individual polyelectrolytes, for ex- 
ample, the effective charge or the bending rigidity. The 
screening also determines forces acting between polyelec- 
trolytes and both thermodynamic and transport proper- 
ties of their solutions. Here we concentrate on a rigid 
polyelectrolyte with a fixed charge distribution. Two 
standard problems are considered below — an infinite 
flat surface and an infinite cylinder — each uniformly 
charged with the surface density —a < 0. The stan- 
dard approach for a description of such problems is the 
Poisson-Boltzmann equation (PBE) for the selfconsistent 
electrostatic potential ^>(r) 

**--££**.•*(-$)■ (1) 

Here e is the charge of a proton, D ~ 80 is the dielec- 
tric constant of water, i^e is the charge of a small ion of 
sort i and Nqi is their concentration at the point where 
■0 = 0. The number of papers using the analytical and 
numerical solutions of Eq. (Q) is extremely largea. On 
the other hand there is understanding that Eq. (fy) ne- 
glects ion-ion correlations and is not exact. Deviations 
from the distribution of cfewge predicted by PBE were 
demonstrated numericallyEra for the following problem. 
Consider screening of a charged surface, x — 0, of a mem- 
brane or a film by a water solution occupying halfspace 
x > 0. Assume that there is only one sort of counterions 
with charge Ze > and their concentration N(x) — > 
at x — > oo. In this case solution of Eq. (nl) is very simple 
and has the Gouy-Chapman form 



N(x) 



1 



2nl{X + xy 
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where A = Ze/(2irla) is the Gouy-Chapman length, 
I = Z 2 l B and l B = e 2 /(Dk B T) ~ 0.7 nm is the Bjer- 
rum length. At large Z and <j, the length A can become 
of the order of the size of the water molecule or even 
smaller. For example, at Z = 3 and a = 1.0 e/nm 2 length 
A = 0.08 nm. Thus, almost all counterions are located in 
the first molecular layer at the surface or, in other words, 
they condense at the very surface of polyelectrolyte. This 
raises question about a role of their lateral correlations 
ancLvalidity of Eq. (0) . It was found by numerical meth- 
odsQ that for a typical charge density a deviations from 
Eq. (g) are not large for monovalent counterions, but they 
strongly increase with the charge of counterions Z . It was 
suggested in Ref. [|, % %■, % that at Z > 2 repulsion be- 
tween multivalent counterions condensed at the surface 
is so strong that they form a two-dimensional strongly 
correlated liquid (SCL) in which the short order of coun- 
terions is similar to that of Wigner crystal (WC). This 
idea was used to demonstrate that two charged surfaces 
in the presence of multivalent counterions can attract 
each other at small distances. 

The goal of this paper is to develop a simple analytical 
theory of the influence of SCL on a screening atmosphere 
of a single charged surface. It is shown below that the co- 
hesive energy of SCL provides additional binding of ions 
to the surface. PBE fails to describe this correlation ef- 
fect. At the same time, PBE works well far from the sur- 
face where N(x) is small and correlations are not impor- 
tant. To describe screening at large distances we derive a 
new boundary condition for N(x) at x = 0, which takes 
the effect of SCL into account. Then we solve PBE with 
this boundary condition for standard problems of screen- 
ing of charged flat surface and a cylinder for different 
salt concentration N(oo) in the bulk of solution. In the 
case of a cylinder, we show that in practical conditions 
the conventional picture of the Onsager-Manningu con- 
densation should be strongly modified when dealing with 
multivalent ions. Since the counterions are tighly bound 
to the surface of cylinder the net linear charge density of a 
screened cylinder is smaller than in the Onsager-Manning 
theory and depends on a bare density. The net charge 



density can also cross zero and become positive. In the 
case of a flat surface similar phenomena are found which 
lead to the screening atmosphere different from Eq. (0) 
and are in qualitative agreement with numerical resultsu. 
Because of simplicity and the use of the universal bound- 
ary condition our theory complements direct numerical 
methodsu which until now were not able to study the case 
of Z > 3 and other than planar geometries. 

Let us present arguments for the existence of SCL and 
then derive a new boundary condition for N(x). As 
we mentioned above and will show below for a typical 
a > le nm~ 2 and Z > 3, almost all charge of the plane 
is compensated by SCL of counterions at its surface, 
which has a two-dimensional concentration almost equal 
to n = a/Ze. The minimum of Coulomb energy of coun- 
terion repulsion and their attraction to the background 
is provided by triangular close packed WC of counteri- 
onst3. At T — the energy per ion of such WC, e(n), can 
be estimated as the energy of attraction, —Z 2 e 2 /DR, of 
an counterion to its Wigner-Seitz cell (correlation hole) 
which is approximately a disc with radius R — (nn)^ 1 / 2 
and charge — Ze. More accurately 



e(n) 



-an 



1 / 2 Z 2 e 2 D~ 



(3) 



where a = 1.96. The inverse dimensionless temperature 
of SCL is usually written in units 



2 D 2 



Z 2 e 



RDk B T 



0.9 



\e(n)\ 
k R T' 
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For example, at a = 1.0 e/nm 2 and room temperature, 
Eq. (|) gives T = 3.5, 6.4, 9.9 at Z = 2, 3, 4. Thus 
for multivalent counterions we deal with low tempera- 
ture situation. T is the large parameter of our theory. In 
its terms R/X ~ 2r > 1 and l/R ~ T > 1. For example, 
at Z = 3 and u = 1.0 e/nm 2 lengths A, R and I are equal 
to 0.08, 1.0, 6jlnm respectively. 

It is knownllil , however, that WC melts at very low 
temperature near T ~ 130. So in the range of our 
interest, 3 < T < 15, we deal with SCL. Thermody- 
namic properties of such SCX*,or one component plasma 
were studied numerically E30. In the large range 0.5 < 
r < 50 excessive internal energy of SCL per counterion, 
e(n,T) = k B Tf(T), was fitted by the expressior£3 



f(T) = -l.lF + 0.58r i/4 - 0.26 



(5) 



with an error less than 2%. The first term on the right 
side of Eq. (||) dominates at large T and leads to Eq. (||). 
The other two terms provide a relatively small correction 
to the energy of WC. It is equal to 11% at T = 5 and 
to 5% at T = 15. The reason for a such small differ- 
ence is that short range order in SCL is similar to that of 
WC. For the free energy per unit area, F, we can write 
F = F{T = 0.5) + nk B T f[ 5 f(T')dr'/T', so that for the 
chemical potential which we need below to describe the 
equilibrium of SCL with the gas-like phase we obtain 



(J.(n,T) = -k B Tln(n w /n) + u s +/i c (n,T), (6) 

Mc (n, T) = -k B T(lMT - 2.6ir 1/4 + 0.26 lnT + 0.13). (7) 

Here /i c is contribution of correlations to the chemical 
potential and we replaced n(T = 0.5) by the chemical po- 
tential —k B Tln(n w /n) + fj, s of an ideal two-dimensional 
solution of ions in the surface layer of water with two- 
dimensional concentration n w . The term u s is the hy- 
dration free energy per ion at the surface and at n <C n w 
does not depend on the concentration of ions n. The first 
term of Eq. (M) corresponds to zero temperature WC and 
can be found directly from Eq. (ph. 

We show below that when an counterion moves away 
from SCL, it leaves behind its negatively charged corre- 
lation hole. We will also see that the potential energy 
of attraction to this hole becomes smaller than k B T at 
x > l/A. On the other hand the selfconsistent potential 
of strongly screened surface is so small that it changes 
by k B T only at exponentially large length A, which is 
defined below. Therefore the condition of equilibrium 
between SCL and the layer 1/4 <C x ^C A is 



n(n) = u(N), 
where u(n) is given by Eq. (g) and 

H(N) = -k B T\nN w /N + m, 



(8) 



(9) 



is the chemical potential of the bulk gas-like phase, N w 
is the bulk concentration of water and /if, is the bulk hy- 
dration free energy which does not depend on N. Using 
Eq. (||) and Eq. (||) and solving Eq. (||) for N we obtain 
at l/A < x < A 



N{Q) = - exp 

W \ KrJ 



(10) 



where w = (n w /N w ) exp[(fii, — fi s )/k B T)]. Below we as- 
sume for simplicity that /if, = u s . In this case w is the 
length of the order of size of the water molecule (for es- 
timates we use w = 0.3 nm). 

The notation N (0) reflects the fact that this value plays 
the role of a new boundary condition at i C A for so- 
lution of PBE which as we will see has large character- 
istic length A. Due to the large value of \u c (n,T)\, the 
concentration A^(0) can be very small. For example, at 
Z = 3 and a = 1.0 e/nm -2 , when T — 6.4 and according 
to Eq. (0) \u c (n,T)\/k B T = 7.0 we obtain N(0) = 10 24 
m -3 = 1.6 mM. Switching to Z = A we get V = 9.9, 
\Ha(n,T)\/k B T = 12.4 and N(0) = 5.5 /xM. It is clear 
now that the \fx c {n, T)\ plays in our problem the role 
similar to the work function for thermal emission or to 
the free energy of chemisorption. Thus we see that corre- 
lation effects in SCL provide nonspecific strong binding 
of counterions to the macroion surface. Qualitative ar- 
guments for such binding can be found in Ref. H 

We would like to stress that such binding does not 
happen at Z = 1. Indeed, at a = 1.0 e/nm -2 one 



obtains from Eq. (g) and Eq. (0) that T = 1.2 and 
fi c (n,T)/k B T — 0.4. Therefore the boundary condition 
Eq. (|1C]) does not produce nontrivial effects and standard 
solutions of PBE remain approximately valid. 

Below we justify the role of the distance 1/4 and 
give an idea how N(x) evolves from n/X at x ~ A to 
N(0) at x = 1/4. Let us move one ion of the crys- 
tal along the x axis. As we mentioned above, it leaves 
behind its correlation hole. In the range of distances 
A <C x *C R, the correlation hole is approximately the 
disc of charge with radius R (the Wigner-Seitz cell) and 
the ion is attracted to the surface by its uniform electric 
field E = 2tt(t/D. Therefore, if A were larger than w we 
would get N(x) — (n/X) exp(—x/\) at x <C R. In the 
cases of our interest X < w and N = n/w at x < w, 
while at w <C x <C R 



N(x) = — exp(-x/A). 

w 



(11) 



At x 3> R the correlation hole radius grows and becomes 
of the order of x. Indeed SCL can be considered as a 
good conductor in the plane (y,z). It is known that a 
charge at distance x from a metallic plane attracts an 
opposite charge into a disc with the radius ~ x or, in 
other words, creates its point like image on the other side 
of the plane at the distance 2x from the original charge. 
The same thing happens to SCL. The removed ion repels 
other ions of SCL and creates a correlation hole in the 
form of the negative disc with the charge —Ze and the 
radius x. The correlation hole attracts the removed ion 
and decreases its potential energy by the Coulomb term 
U(x) = —Z 2 e 2 /4Dx. This effect provides the correction 
to the activation energy of N(x) : 



N(x) — — exp 
w 



\fj,(n)\ ~Z 2 e 2 /4Dx 
kWT 



(x>R). (12) 



The similar "image" correction to the work function of a 
metal is well-known in the theory of thermal emission. 
The correction decreases with x and at x = 1/4, becomes 
equal to fcgT, so that N(x) saturates at the value N(0). 
The dramatic difference between the exponential decay 
of Eqs. @, (g||) and 1/x 2 law of Eq. @ is obviously re- 
lated to correlation effects neglected in PBE. Recall that 
it was assumed in the beginning of this paper that the 
charge of the surface is almost totally compensated by 
SCL. Exponential decay of N(x) with x confirms this as- 
sumption and at P 3> 1 makes this theory self-consistent. 
Consider now what happens with N(x) at distances 
x ^> 1/4. At such distances, correlations of the removed 
ion with its correlation hole in SCL are not important 
and correlation between ions of the gas phase are even 
weaker because N(x) is exponentially small. Therefore, 
one can return to PBE. Solution of PBE with the bound- 
ary condition ( |l0| ) and JV(oo) = is similar to Eq. (||): 



N{x) 



2nl (A + x) 2 



(x » 1/4), 



(13) 



where the new renormalized Gouy-Chapman length, A, 
is exponentially large 



A 



< 2 " w <°»" /2 - i/ss 



exp 



W(n,T)\ 
2k B T 



(14) 



For example, at a = 1.0 e/nm -2 Eq. (ph gives A ~ 
0.8,5.3,68 nm at Z — 2,3,4. Comparing these lengths 
with 1/4 = 0.7,1.6,2.8 nm respectively we see that 
A > 1/4 for Z > 3, what justifies the use of Eq. @ 
as the boundary condition for the large distance solution 
of PBE. At Z = 2, however, A ~ 1/4 and our theory 
works only qualitatively. This is the reason why we do 
not compare our N(x) quantitatively with numerical re- 
sultsa, obtained only for Z = 2. 

Using Eq. (|l^) one finds that the total surface charge 
density located at distances x < 1/4 



- v / N(0)/(2ttI b ) = -<r(A/A). 



(15) 



The exponentially small a* is a result of correlation ef- 
fects, which strongly bind countcrions to the surface. 
Corrections to fi c (n,T) and N(0) related to such small 
a* can be neglected. Until now we talked about the case 
when N(oo) = 0. Let us now assume that there is a con- 
centration N(oo) of Z:\ salt in the bulk of solution, so 
that in the bulk N(oo) of counterions is neutralized by 
N-(oo) = ZN(oo) of monovalent coions. This adds the 
Debye-Hiickel screening radius 



r s = (4TrlN(oo)(l + Z- 2 )) 



-1/2 



(16) 



to the problem. If N(oo) -C iV(0) screening radius 
r s ^> A and the fact that N(oo) is finite changes only 
the very tail of Eq. ( |l3| ) making decay of N(x) at x ^> r s 
exponential. At i « r s still N(x) ^> ZN-(x) and all 
previous results are valid. However, when iV(oo) ap- 
proaches A(0), the solution changes dramatically and 
a* vanishes. Indeed when N(oo) = N(0) concentra- 
tion N(x) = N(oo)cxp(-Zeip/k B T) = N(0) stays con- 
stant and potential ip(x) — at x > 1/4. This means 
that the surface is completely neutralized at distances 
< x < 1/4. 

If A(oo) ^> A(0) negative charges dominate at 
x <C r s . Indeed in the PBE approach, N(x) = 
N(oo)exp(-Zei}/k B T),N-(x) = A_(oo) exp(eip/k B T) 
and when concentration N(x) decreases iV_(x) increases. 
One can derive a boundary condition for N—(x) at x = 
from these equations 



M(0) = ZN(oo)[N(oo)/N(0)} 1/z , 



(17) 



where N(0) is given by Eq. ([io|). Then the solution 
of PBE for N-(x) at a; C r s has the form simi- 
lar to Eq. (JH) N_(x) = (27t/b) _1 ( A - + x )~ 2 > where 
A_ = (27r/ B A_(0)) -1 / 2 . To compensate the bulk neg- 
ative charge the positive surface charge density of SCL 
becomes larger than a, so that the net surface charge 
density a* > 0. Similarly to Eq. (Ill), it is 



a* =e y/N_(0)2Trl B = e/2irl B A-- (18) 

This phenomenon is called charge inversion and is, of 
course, impossible in the framework of the standard PBE. 
Technically, charge inversion follows from the small value 
of iV(0) in Eq. ([10|). Its physics is related to the strong 
binding of counterions at the charged surface due to for- 
mation of SCL. Remarkably, when T 3> 1, this phe- 
nomenon happens under the influence of very small con- 
centration of salt. 

Let us switch to the screening of an infinite rigid cylin- 
der with a radius a, a surface charge density — a or, in 
other words, with a linear charge density rj = —2'Kaa. 
We assume that a is large enough so that the surface of 
the cylinder is covered by a two-dimensional SCL with 
R < 2ira and with F>1. (This means that the follow- 
ing important inequality, \rj\ 3> r\ c = eZ/l, is satisfied). 
Such a cylinder can be a first order approximation for 
the double helix DNA, where a = 1 nm, r\ — 5.9 e/nm, 
a = 0.94 e/nm 2 , and at Z = 3 length R ~ 1 nm and 
I — 6.3 nm. Below we assume that, as in this example, 
I » a. 

A screening atmosphere of a cylinder is described by 
concentration N(r), where r is the distance from the 
cylinder .axis. For \rj\ 3> r\ c the solution of PBE is 
knownl!3' : to confirm the main features of the Onsager- 
ManningS picture of the counterion condensation. The 
screening charge \r\\ — rj c is located at the cylinder surface, 
while the rest of the screening charge, rj c , at JV(oo) = 0, 
is spread in the bulk of the solution. This means that 
the net charge density of the cylinder, rf , equals — rj c 
and does not depend on rj. At a finite N(oo) the charge 
density rf is screened only at linear screening radius r s . 

Does SCL at the surface of the cylinder change these 
elegant statements? We show below that 

ln[Y(0)/Y(oo)] 
1 " /c ln[4/(7r7V(oo)Z 3 )]' V ; 

It is clear from Eq. ( |l9| ) that if two logarithms are close 
to each other, i. e. if ln(N (0) 2 1 3 / N (oo)) > 1 the 
Onsager-Manning theory is approximately correct and 
rf approaches — rj c . Note, however, that concentration 
N(0) itself is exponentially small so that values of iV(oo) 
at which rf is close to —rj c are unrealistically small. 
On the other hand, in disagreement with the Onsager- 
Manning theory one obtains from Eq. ( |l9|) that \rf\ <C r\ c 
when concentration N(oo) of the salt is still exponentially 
small, namely N(0) 2 l 3 < iV(oo) < Y(0). Moreover, ac- 
cording to Eq. ( p9J ) rj* vanishes at N(oo) = N(0). This 
result is easy to understand without calculations. In- 
deed, in this case N(r) = iV(oo) exp(— Zeijj(r) / ksT) = 
iV(0) stays constant and ip( r ) = at all r > 1/4, so 
that all charge of the polyelectrolyte is compensated in- 
side cylinder with r = 1/4. The difference from the 
Onsager-Manning theory becomes even more apparent 
at N(oo) > N(0) when the density rj* becomes posi- 
tive, once more demonstrating the charge inversion. Note 
that the charge inversion takes place still at exponentially 
small N(oo). One can show that a positive rj* continues 



to grow with growing iV(oo) until it reaches critical den- 
sity e/ls = rj c / ' Z at which Onsager-Manning condensa- 
tion of monovalent coions begins. Charge rj* includes all 
counterions for which the binding energy exceeds ksT. 
Therefore, anomalous charge density of a polyelectrolyte 
•q* can be measured in electrophoresis experiment. Fi- 
nally, Eq. ([l9|) shows that rj* does depend on bare charge 
density rj trough n in Eq. (10). Thus, for screening by 
multivalent ions at T ;§> 1 and at reasonably large iV(oo) 
all statements of Rcf. g, [F| are qualitatively incorrect. 

To derive Eq. (|l9|), we focus at distances 1/4 < r < r s , 
where the electrostatic potential of the linear charge den- 
sity if is not screened and the boundary condition of 
Eq. ( |1C)| ) can be used to write 



~ N(0) exp (^ ln(4r/0 



(20) 



At r = r s concentration N(r s 
equation for rj* we get 



iV(oo). Solving this 



V = -Vc 



ln [AT(0)/iV(oo)] 
ln(4r s /0 2 



(21) 



Finally, using Eq. (jlq), we arrive at Eq. ( |l9|) 

We are grateful to V. A. Bloomfield, M. M. Fogler, A. 
Yu. Grosberg, R. Kjellander, S. Marcelja, R. Podgornik 
and I. Rouzina for valuable discussions. This work was 
supported by RFSF (V. P.) and by NSF DMR-9616880. 



1 Permanent address: A. F. Ioffe Physico- Technical Institute, 
194021, St. Petersburg, Russia. 

see, for example, M. D. Frank-Kamenetskii, V. V. Ansele- 
vich, and A. V. Lukashin, Sov. Phys. Usp. 30, 317 (1987). 

3 L. G. Gulbrand, Bo Jonsson, H. Wennerstrom, and P. 
Linse, J. Chem. Phys. 80, 2221(1984). 

4 R. Kjellander and S. Marcelja, Chem. Phys. Lett. 114, 
124(E) (1985); R. Kjellander, Ber. Bunsenges. Phys. Chem. 
100, 894 (1996); 

5 I. Rouzina and V. A. Bloomfield, J. Phys. Chem. 100, 9977, 
(1996). 

6 N. Gronbech- Jensen, R. J. Mashl, R. F. Bruinsma, and W. 
M. Gelbart, Phys. Rev. Lett. 78, 2477 (1997) 



7 J. J. Arenzon, J. F. Stilck, and Y. Levin, cond 



mat/9806358 . 
s B. I. Shklovskii, Phys. Rev. Lett. 82, 3272 (1999). 
9 G. S. Manning, J. Chem. Phys. 51, 924 (1969). 
10 L. Bonsall, A. A. Maradudin, Phys. Rev. B15, 1959 (1977) 

R. C. Gann, S. Chakravarty, and G V. Chester, Phys. Rev. 

B 20, 326 (1979). 

12 H. Totsuji, Phys. Rev. A 17, 399 (1978). 

13 B. Zimm, M. LeBret J. Biomol. Struct. Dyn. 1, 461 (1983). 



